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1.  Summary  of  Research  Progress 

This  section  summarizes  progress  made  during  the  reporting  period  in 
two  major  areas  of  the  research.  It  begins  with  a  brief  introduction  and 
overview  of  the  research  program  in  §1.1.  Following  this,  §1.2  gives  a  summary 
of  results  from  our  study  of  the  non-equilibrium  structure  of  combustion  in 
turbulent  flows.  Next,  §1.3  summarizes  results  from  measurements  of  the  large 
scale  topology  and  scaling  properties  of  the  scalar  dissipation  rate  layers  which 
form  the  basic  feature  of  the  molecular  mixing  process,  and  shows  results  from 
comparisons  with  an  analytical  model  of  the  physical  process  that  leads  to  this 
structure. 

1.1  Introduction  and  Overview 

This  research  program  in  the  Gas  Dynamics  Laboratories  at  The 
University  of  Michigan  consists  of  a  combined  experimental  and  theoretical 
investigation  into  the  structure  of  combustion  reactions  in  turbulent  flows,  with 
the  following  three  major  objectives: 

(i)  to  develop  and  implement  rtew,  high-resolution,  multi-dimensional, 
quantitative,  imaging  capabilities  for  obtaining  direct  experimental 
measurements  of  the  fine  scale  structure  associated  with  the 
molecular  mixing  and  chemical  reaction  processes  in  turbulent  flows, 

(ii)  to  use  these  new  measurement  techniques  to  experimentally 
investigate  the  essential  physical  characteristics  of  the  fine  stmcture 
of  molecular  mixing  and  chemical  reactions  in  turbulent  flows,  and 

(Hi)  to  incorporate  results  from  these  new  experimental  measurements 
into  an  improved  understanding  of  the  molecular  mixing,  chemical 
reaction,  and  local  extinction  processes  in  reacting  turbulent  flows  of 
interest  for  practical  air-breathing  propulsion  systems. 

The  overall  goal  of  this  program  is  to  contribute  to  the  development  of 
new  models  for  turbulent  combustion  in  air-breathing  propulsion  systems. 
Such  models  are  viewed  as  an  essential  tool  in  the  development  of  future 
propulsion  systems.  Their  development  is,  in  turn,  limited  by  the  current  lack 
of  understanding  of  the  nature  of  molecular  mixing  in  turbulent  flows  and  its 
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relation  to  complex  combustion  chemistry.  This  program  is  yielding  important 
insights  into  the  structure  of  reacting  turbulent  flows  at  the  molecular  diffusion 
and  chemical  reaction  scales  of  the  flow.  It  is  based  on  direct  measurements  of 
the  fully  space-  and  time- varying  fine  structure  of  conserved  scalar  fields  in 
nonreacting  turbulent  flows,  and  the  relation  between  these  and  the  physical 
structure  of  the  corresponding  chemical  species  concentration  fields  K,(x,r)and 
the  associated  chemical  reaction  rate  fields  in  nonequilibrium  reacting 

turbulent  flows. 

This  Annual  Report  focuses  on  results  obtained  during  the  first  year  of 
this  Phase  II  research  effort.  A  number  of  basic  concepts  essential  to  the 
approach  taken  here,  as  well  as  results  from  our  initial  two-year  start-up  effort, 
are  contained  in  the  1989-91  Final  Report  (see  also  AFOSR-TR-91-0111).  Results 
have  also  been  described  in  greater  detail  in  several  technical  articles  appearing 
in  the  refereed  archival  literature*  as  well  as  in  several  monograph  articles  and 
conference  presentations*. 

Work  in  this  program  is  proceeding  along  two  major  avenues.  The  first 
involves  measurements  conducted  at  moderate  Reynolds  numbers  and  large 
Schmidt  numbers  using  a  unique,  four-dimensional,  high-resolution,  non- 
intrusive  imaging  capability  developed  under  our  Phase  I  effort.  These  four¬ 
dimensional  imaging  measurements  are  yielding  fully-resolved,  giga-byte  sized, 
experimental  data  volumes  of  the  precise  conserved  scalar  field  structure 


'  see  Dahm,  Southerland  &  Buch  (1991)  Physics  of  Fluids  A  3,  1115-1127. 

Tryggvason  &  Dahm  (1990)  Combustion  &  Flame  83,  207-220. 

Chang,  Dahm  &  Tryggvason  (1991)  Physics  of  Fluids  A  3,  1300-1311. 

Southerland,  Porter,  Dahm  &  Buch  (1991)  Physics  of  Fluids  A  3,  1385-1392. 

Dahm,  Su  &  Southerland  (1992)  Physics  of  Fluids  A  4,  2191-2206. 

*  see  Dahm,  Southerland  &  Buch  (1991)  in  Applications  of  Laser  Technigue.s  to  Fluid 
Mechanics.  R.  Adrian,  Ed.,  Springer  Verlag,  Berlin. 

Dahm  &  Buch  (1990  in  Chemical  Reactions  and  Physical  Processes  in  Turbulent  Flows. 
J.  Hunt,  Ed.,  Cambridge  University  Press. 

Buch,  Dahm,  Dibble,  &  Barlow  (1992)  in  Proceedings  of  the  24th  International 
Symposium  on  Combustion,  The  Combustion  Institute,  Pittsburgh. 

Dahm,  W.J.A  (1992)  in  Proceedings  of  the  13th  Symposium  on  Turbulence.  University  of 
Missouri-Rolla,  Rolla,  MO. 

Dahm,  W.J.A  (1992)  in  Turbulence  and  Molecular  Processes  in  Combustion.  T.  Takeno, 
Ed.,  Elsevier. 
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on  the  inner  scales  of  turbulent  flows  (see  Fig.  1).  These  data  are  simultaneously 
resolved  in  all  three  space  dimensions  and  in  time,  allowing  direct  extraction  of 
the  instantaneous  scalar  energy  dissipation  rate  field  V^-V^(x,r)  and  its  time 
evolution.  These  two  fields  can  be  used  to  deduce  the  mixing  and  chemical 
reaction  processes  in  turbulent  reacting  flows.  Data  of  this  type  have  previously 
been  conceivable  only  from  large-scale  direct  numerical  simulations  (DNS)  of 
the  full  Navier-Slokes  equations,  yet  these  experimental  measurements  are  now 
possible  in  real  turbulent  flows  and  can  reach  conditions  outside  the  range  of 
accessibility  of  such  numerical  studies.  Under  our  Phase  I  effort,  these 
measurements  were  used  to  show  that  essentially  all  of  the  molecular  mixing  in 
turbulent  flows  occurs  in  thin,  strained,  laminar  diffusion  layers  which  can  be 
accurately  modeled  by  a  simple  canonical  solution  of  the  Navier-Stokes  and 
conserved  scalar  transport  equations.  Under  the  present  Phase  II  effort,  the 
larger  scale  topology  and  scaling  properties  of  these  dissipation  layers  are  being 
investigated  to  develop  a  complete  model  of  the  mixing  process.  These 
processes  are  driven  by  the  continual  stretching  and  folding  imposed  on  the 
dissipation  layers  by  the  velocity  gradient  field  in  the  underlying  turbulent  flow. 
Scaling  properties  deduced  thus  far  show  a  remarkable  similarity  in  the  resulting 
dissipation  field  structure  between  high  Reynolds  number  turbulent  flows  and 
simple  low  Reynolds  number  chaotic  flows.  This  suggests  great  potential  for  the 
development  of  simple  models  of  the  molecular  mixing  process  in  practical 
turbulent  flows. 

Complementing  these  studies  of  molecular  mixing  in  turbulent  flows,  a 
second  major  effort  of  our  work  involves  highly-resolved,  multi-dimensional 
imaging  measurements  of  5c  «  1  conserved  scalar  mixing  in  gaseous  turbulent 
flows.  As  is  the  case  in  our  large  Sc  measurements,  the  highly  resolved  nature  of 
the  measurements,  together  with  the  very  high  signal  quality  achieved,  permits 
direct  differentiation  of  the  data  to  allow  accurate  evaluation  of  derivatives 
involved  in  the  scalar  gradient  field  V^(x,r).  This  yields  the  structure  of  scalar 
energy  dissipation  rate  field,  permitting  highly  detailed  analysis  of  the  structure 
of  chemical  reactions  in  nonequilibrium  combustion  processes  occurring  in 
turbulent  flows.  The  local  instantaneous  thermochemical  state  of  the  reacting 
turbulent  flow  can  be  specified  from  a  two-parameter  nonequilibrium 
formulation  in  tein«>  of  the  local  conserved  scalar  and  scalar  dissipation  rate 
values.  This  procedure  is  used  to  examine  the  reaction  zone  structure  in 
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turbulent  flames  under  varying  degrees  of  chemical  nonequilibrium.  These 
measurements  are  providing  a  previously  inaccessible  level  of  detailed 
information  about  the  structure  of  combustion  processes  in  turbulent  flows. 
Results  obtained  permit  the  first  direct  view  into  the  structure  of  the  species 
chemical  reaction  rate  fields  in  reactive  turbulent  flows  under  various 

degrees  of  chemical  nonequilibrium. 

1.2.  Non-Equilibrium  Structure  of  Turbulent  Combustion 

Experimental  studies  to  identify  the  canonical  structure  of  the 
instantaneous  chemical  reaction  rate  fields  Wj(x,r)  in  turbulent  diffusion  flames 
play  a  key  role  in  the  development  of  physically-based  models  of  combustion  in 
turbulent  flames.  Combustion  models  for  turbulent  reacting  flow  ultimately  aim 
to  represent  the  connection  between  the  chemical  source  terms  on  the  right 
side  of  the  species  conservation  equations  and  the  turbulent  dynamics  of  the 
underlying  flow  field.  Studies  of  the  reaction  zone  structure  under  various 
conditions  of  turbulent  flow  are  typically  centered  around  single-point  or  planar 
imaging  measurements  of  certain  species  concentration  fields  Tjfx.r),  from 
which  the  reaction  rate  field  stmcture  is  inferred.  Instantaneous  concentrations 
of  such  species  as  OH  or  CH,  for  example,  have  often  been  interpreted  as 
approximately  marking  the  reaction  rate  field  stmcture.  These  concentration 
fields,  however,  reflect  the  cumulative  effects  of  turbulent  transport,  molecular 
diffusion,  and  chemical  reaction,  and  thus  can  give  a  potentially  misleading 
picture  of  the  underlying  chemical  reaction  rate  field  stmcture 

Direct  measurements  of  the  reaction  rate  fields  are  not  possible,  but  they 
can  be  deduced  from  instantaneous  measurements  of  a  conserved  scalar  field 
^(x,r).  In  the  limit  of  chemical  equilibrium,  for  example,  the  stmcture  of 
individual  species  reaction  rate  fields,  w,(x,r),  are  relatively  simply  related  (Bilger 
1976)  to  the  instantaneous  conserved  scalar  value  and  the  instantaneous 
dissipation  rate,  (^gScr'V^-V^fx.r),  of  the  scalar  energy  per  unit  mass, 

Away  from  equilibrium,  the  connection  between  the  scalar  field  stmcture  and 
the  reaction  rate  fields  is  not  as  simple.  One  approach  to  this  is  to  use  a  strained 
laminar  diffusion  and  reaction  layer  formulation  to  represent  the  interaction 
between  the  turbulent  flow  dynamics  and  the  nonequilibrium  reaction 
chemistry. 
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We  are  using  this  approach  to  examine  the  structure  of  instantaneous 
species  mass  fraction  fields  and  the  associated  instantaneous  chemical  reaction 
rate  fields  in  a  turbulent  jet  diffusion  flame  under  '.aiying  degrees  of 
nonequilibrium  chemistry.  The  objectives  are,  firstly,  to  assess  the  underlying 
canonical  structure  of  nonequilibrium  reaction  rate  fields  in  turbulent  diffusion 
flames  and,  secondly,  to  determine  the  extent  to  which  certain  species 
concentration  imaging  measurements  can  be  used  as  indirect  indicators  of  the 
reaction  field  structure.  The  presenr  results  are  for  OH  and  NO  mass  fraction 
and  reaction  rate  fields  for  hydrogen-air  chemistry.  The  and  Wj(x,r)  fields 

are  constructed  from  fully-resolved  imaging  measurements  of  the  instantaneous 
conserved  scalar  field  ^(x,r),  obtained  via  planar  Rayleigh  scattering  from  the 
self-similar  far  field  of  a  nonreacting  turbulent  jet  flow. 

Details  of  the  flow  facility,  operating  conditions,  and  imaging  diagnostics 
are  described  elsewhere^  For  the  present  purposes,  the  key  points  are  that,  at 
the  flow  conditions  selected  and  for  the  imaging  set-up  used,  the  signal-to-rms- 
noise  ratio  that  exceeded  100:1,  and  the  smallest  conserved  scalar  gradient  scales 
were  fully  resolved  in  both  space  and  time,  reaching  below  the  finest  local  scalar 
gradient  lengthscale  in  the  turbulent  scalar  field.  This,  together  with  the  signal 
quality  achieved,  allows  differentiation  of  the  measured  scalar  field  ^(x,r)  to 
produce  the  scalar  dissipation  rate  field  V^-V^(x,r).  The  self-similar  far-field 
scaling  of  this  flow  allows  these  scalar  and  dissipation  rate  data  to  be  mapped  to 
any  downstream  location.  The  data  are  then  coupled  with  a  strained  laminar 
flamelet  library  for  hydrogen-air  chemistry  to  produce  the  instantaneous 
species  concentration  and  reaction  rate  fields  for  several  increasing  levels  of 
chemical  nonequilibrium.  The  Reynolds  number  is  held  constant  while  the 
Damkohler  number  is  systematically  varied  by  changing  the  effective  nozzle 
diameter.  Comparisons  at  different  downstream  locations  in  the  jet  allow  effects 
due  to  the  decreasing  scalar  values  to  be  separated  from  effects  due  to  the 
varying  degree  of  nonequilibrium  chemistry.  We  apply  this  approach  to 
examine  the  OH  and  NO  mass  fraction  and  reaction  rate  fields  under  conditions 
reaching  to  deep  nonequilibrium. 


*  see:  Buch,  Dahm,  Dibble  &  Barlow  (1991)  Structure  of  equilibrium  reaction  rate  fields  in 
turbulent  jjet  diffusion  flames.  Proceedings  of  the  24th  (International)  Symposium  on 
Combustion,  Sydney,  Australia,  July  5-10,  1992.  See  also  Buch  &  Dahm  (1991)  Fine  scale 
structure  of  conserved  scalar  mixing  in  turbulent  shear  flows.  Sc  »  1,  Sc  =  1,  and  implications 
for  reacting  flows.  University  of  Michigan  Report  No.  026779-5- 
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In  our  previous  year  report  we  showed  a  typical  256  x  512  measured  data 
plane  of  the  instantaneous  conserved  scalar  field  structure  ^(x,r).  The  256 
different  colors  denoted  ranges  of  conserved  scalar  values,  with  pure  blue 
beginning  at  ^(x.r)  -  0,  corresponding  to  pure  air,  and  increasing  uniformly  to 
pure  red  denoting  the  highest  conserved  scalar  values  in  the  data.  The  jet 
centerline  ran  down  the  right  edge  of  the  data  plane,  and  the  axes  indicated  the 
spatial  extent  of  the  data  in  terms  of  the  local  strain-limited  scalar  gradient 
lengthscale  Xq.  Derivatives  of  the  measured  conserved  scalar  field  values  give  the 
projection  of  the  true  three-dimensional  scalar  gradient  vector  ViJ(x,r)  into  the 
plane.  We  showed  the  logarithm  of  the  corresponding  scalar  dissipation  rate 
field  formed  from  this  projection,  The  256  different  colors 

denoted  the  local  scalar  dissipation  rate,  with  black  beginning  at  -  0,  and 

pure  blue  through  pure  red  denoted  logarithmically  increasing  values  as 
indicated.  These  figures  should  be  referred  to  in  the  following  discussion. 
While  two-dimensional  dissipation  values  underestimate  the  magnitude  of  the 
full  dissipation  rate  wherever  the  scalar  gradient  vector  has  a  strong  out-of-plane 
component,  the  present  values  nevertheless  provide  an  accurate  picture  of  the 
topological  structure  of  the  dissipation  field. 

The  local  outer  variables  «(x)  and  5(x),  and  the  mass  fraction-based 
conserved  scalar  field  in  the  far  field  of  axisymmetric  turbulent  jets  follow 

known  self-similar  scalings  with  downstream  location  x.  This  allows 
measurements  such  as  ours  to  be  uniformly  rescaled  to  map  the  instantaneous 
scalar  and  dissipation  rate  fields  to  any  x-location  in  the  far  field  of  the  flow.  A 
set  of  hypothetical  chemical  reactions  occurring  in  this  conserved  scalar  field 
would  then  produce  a  turbulent  jet  diffusion  flame  of  length  L,  where  the  flame 
tip  is  taken  to  occur  at  the  x-location  for  which  the  maximum  conserved  scalar 
value  achieves  stoichiometry  for  the  fuel  and  oxidizer  combination  being 
considered.  The  mass  fraction  and  reaction  rate  fields  for  any  chemical  species  i 
in  this  flame  can  then  be  constructed  from  the  local  values  of  the  scalar  and 
dissipation  rate  fields  using  a  steady  flamelet  approximation. 

Here  we  use  a  flamelet  library  for  pure  hydrogen-air  strained  laminar 
diffusion  and  reaction  layers,  generated  from  numerical  calculations  (supplied  by 
J.-Y.  Chen)  using  a  subset  of  the  Miller  &  Bowman  (1989)  hydrogen-air  kinetic 
mechanism.  Constructing  the  library  requires  mapping  each  (^,  Vi^  VQ  pair  to 
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the  corresponding  strain  rate  e.  The  mapping  obtained  from  the  numerical 
calculations  is  shown  in  Fig.  2a.  For  small  values  of  ^  ,  the  calculations  produce  a 
nonunique  mapping,  presumably  due  to  a  nonuniqueness  in  defining  the 
conserved  scalar  that  results  from  differential  species  diffusivities.  [For  the 
highest  strain  rate  (e  *■  14,000  s~*),  the  map  was  also  nonunique  over  the  entire 
range  of  scalar  values.)  For  this  reason,  the  mapping  was  instead  based  on  the 
analytical  result  for  a  nonreacting  layer,  given  by 

2 

■no  . 

where  f(Q  is  chosen  to  fit  the  numerical  result,  accounting  for  effects  of  heat 
release  and  temperature  dependent  transport  properties.  This  yields  the 
mapping  shown  in  Fig.  2b. 

Here  0  and  0  represent  the  scalar  endpoint  values  on  eithe--  side  of  the 
layer.  For  a  layer  separating  pure  fuel  and  air,  as  in  the  present  library,  these 
values  are  0  and  1.  However,  as  Fig.  1  demonstrates,  none  of  the  layers  in  the  jet 
far  field  are  bounded  by  pure  fuel  and  air,  and  instead  generally  involve  various 
degrees  of  partial  premixing  on  one  or  both  sides  of  the  layer.  A  detailed 
treatment  could  take  this  into  account  but,  as  is  generally  done  in  flamelet 
models,  in  the  present  study  any  effect  of  varying  endpoint  values  is  ignored. 

The  OH  mass  fraction  library  £)  is  shown  in  Fig.  3a,  and  the  reaction 
rate  library  e)  in  Fig.  3b.  A  similar  library  is  constructed  for  NO.  For 

strain  rates  above  e  =  14,000  s~^  the  layer  extinguishes  and  the  flamelet  library 
for  e)  in  Fig.  3b  appropriately  contains  zero  values.  The  mass  fraction 
library  in  Fig.  3a  also  contains  zero-OH  concentrations  in  this  domain  of  YiiO  e), 
but  these  values  are  spurious  since  they  correspond  only  to  the  steady  state 
attained  long  after  extinction.  The  species  concentrations  in  the  transient  layer 
after  extinction  follow  an  exponential  decrease  due  to  the  continued  straining 
action  of  the  flow  field.  A  steady  flamelet  formulation  cannot  account  for  this 
transient  (as  well  as  other  effects  of  unsteadiness),  and  therefore  the  zero  OH 
and  NO  concentrations  resulting  from  extinction  in  the  K,(x,r)  fields  presented 
in  the  following  section  only  denote  where  extinction  has  occurred.  The  results 
obtained  for  the  Wjfx.r)  fields,  however,  do  not  suffer  this  shortcoming,  since 


e  =  TtD 


VC -VC 


exp  2\  erf 


i(C-C) 
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they  represent  inherently  instantaneous  reaction  rate  information. 

Results  are  presented  here  for  two  different  hydrogen-air  turbulent  jet 
diffusion  flames,  having  the  same  far  field  Reynolds  number  -  14,000),  but 
with  conditions  in  the  flame  ranging  locally  from  moderate  to  deep 
nonequilibrium  conditions.  The  jet  conditions  for  each  of  the  cases  considered 
are  shown  in  Fig.  4,  where  the  depth  of  nonequilibrium  is  indicated  relative  to 
the  flame  blowout  limif  for  an  unpiloted  hydrogen-air  turbulent  jet  diffusion 
flame.  With  the  characteristic  temperature  approximately  1800  K,  and  v  ~T^-^, 
this  requires  2  cold-flow  Reynolds  number  of  210,000  to  produce  a  hot-flow 
conserved  scalar  field  consistent  with  the  Re^-  14,000  result  in  Fig.  1. 

Our  results  have  shown  that  Case  (1)  in  Fig.  4  reproduces  the  equilibrium 
fields  given  in  Buch  et  al  (1992).  Cases  (2)  and  (3)  are  used  here  to  examine  the 
effects  of  chemical  nonequilibrium.  Case  (2)  produces  a  turbulent  flame  with  a 
moderate  level  of  nonequilibrium.  It  corresponds  to  a  10  mm  exit  diameter  and 
an  exit  velocity  (after  correction  for  compressibility)  of  nearly  300  m/s.  By 
comparison.  Case  (3)  in  Fig.  4  represents  a  turbulent  flame  under  deep 
nonequilibrium  conditions,  with  an  exit  diameter  of  3  nim  and  exit  velocity  of  1 
km/s,  which  can  be  seen  to  be  only  about  slightly  smaller  than  the  flame 
blowout  limit  (Broadwell  et  al,  1984). 

Figure  5  compares  the  OH  mass  fraction  fields  obtained  for  Cases 

2  and  3.  Both  are  presented  with  the  same  color  assignment  scheme,  as 
indicated  on  the  color  bar,  in  which  pure  blue  denotes  zero  OH  concentration, 
and  colors  increasing  from  blue  to  pure  red  denote  linearly  increasing  OH  mass 
fractions,  clipped  at  0.014.  The  fields  are  shown  at  four  different  downstream 
locations  relative  to  the  mean  flame  length  L,  centered  at  x/L  -  0.25,  0.50,  0.75 
and  1.0.  The  size  of  each  panel  has  been  correctly  scaled  and  placed  for  its 
downstream  location  to  give  a  proper  indication  of  the  structure  of  the  flame. 
Notice  that  the  resulting  OH  fields  look  remarkably  similar  to  direct  OH  planar 
imaging  measurements  in  turbulent  jet  flames. 

In  making  comparisons  between  the  two  cases,  as  well  as  among 
different  panels  within  each  case,  it  is  important  to  note  that  the  degree  of 
nonequilibrium  varies  not  only  between  these  cases,  but  also  with  the 
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Nozzle  Diameter,  do  (mm) 


Fig.  4.  Map  of  source  conditions  for  turbulent  hydrogen-air  jet  diffusion 
flames,  showing  depth  of  chemical  nonequilibrium  relative  to  the  unpiloted 
blowout  limit  (u^/do  =  1.67-106  s-i,).  solid  line  denotes  far-field  hot  flow 
Reynolds  number  Re^  =  14,000  corresponding  to  cold  flow  Rcz  =  210,000. 
Dashed  lines  denote  constant  global  Damkdhler  number  (Uo/ofo  =  const.) 
Shown  are  the  conditions  for  three  cases,  giving  (1)  the  equilibrium  limit, 
(2)  moderate  nonequilibrium,  and  (3)  deep  nonequilibrium. 
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Fig.  5.  Comparison  of  OH 
mass  fraction  fields 
obtained  from  the  con¬ 
served  scalar  imaging  data 
via  the  laminar  flamelet  con¬ 
struction  procedure.  Shown 
are  results  for  Cases  2  (left) 
and  3  (right).  Table  gives 
relative  strain  rates  for  each 
case  at  each  downstream  lo¬ 
cation. 
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downstream  position  x  for  any  given  case.  Owing  to  the  similarity  scalings  5  ~  x 
and  «  ~  xr^  for  this  flow,  the  characteristic  strain  rates  decrease  like  x-2.  Thus 
while  the  strain  rates  in  each  panel  for  Case  3  on  the  right  are  roughly  10  times 
higher  than  in  the  same  panel  for  Case  2  on  the  left  (see  Fig.  4),  the  strain  rates 
among  the  four  downstream  locations  vary  even  more.  For  example,  at  x/L  - 
0.25  the  strain  rates  are  16  times  higher  than  the  values  at  the  flame  tip.  Within 
the  two  cases  for  which  results  are  presented,  the  strain  rates  span  a  range  of 
nearly  200,  as  indicated  in  the  Table  on  the  figure,  which  shows  the  local  strain 
rates  (inverse  local  Damkohler  number)  relative  to  the  value  at  x/L  -  0.25  in  Case 
2. 


Differences  evident  in  the  OH  concentration  fields  in  Fig.  5  when 
comparing  panels  at  the  same  downstream  location  for  the  two  cases  shown  are 
due  solely  to  the  differing  degrees  of  nonequilibrium.  However,  changes  in  the 
Foh(*.0  field  among  panels  corresponding  to  differing  x/L  values  for  the  same 
case  reflect  changes  in  both  the  depth  of  nonequilibrium  as  well  as  changes  in 
the  scalar  and  dissipation  rate  values.  For  the  jet  similarity  scalings,  the  scalar 
field  values  decrease  with  increasing  distance  from  the  jet  source  as  C  ~  x-’. 
Coupled  with  the  linear  increase  in  all  length  scales  in  the  jet  with  increasing  x, 
the  resulting  scalar  dissipation  field  values  decrease  with  distance  from  the  jet 
source  like  V^-V^  ~  x''^.  These  changes  in  both  the  and  V^-V^(x,r)  values 
with  x/L  will  reassign  each  pixel  in  these  data  planes  to  a  different  point  in  the  (^, 
V^-VQ  space  in  Fig.  2,  and  thereby  to  a  different  (^,  e)  point  in  the  flamelet 
library  in  Fig.  3-  Thus,  in  examining  the  results  in  Fig.  5,  as  well  as  those  in  Figs. 
6-8,  it  must  be  borne  in  mind  that  comparisons  made  among  different  cases  but 
at  the  same  downstream  location  show  the  effects  of  the  degree  of 
nonequilibrium  only,  while  comparisons  made  at  different  downstream 
locations  within  the  same  case  demonstrate  the  effects  of  both  nonequilibrium 
level  and  changes  in  the  scalar  field  values.  However,  by  arranging  the  results  in 
the  manner  shown,  it  is  possible  to  separate  the  changes  due  solely  to  the 
degree  of  nonequilibrium  from  those  due  to  changes  in  the  conserved  scalar  and 
dissipation  rate  values. 

With  this  we  can  now  meaningfully  understand  the  nonequilibrium 
structure  of  the  OH  mass  fraction  fields  in  Fig.  5.  Note  firstly  that,  within  Case  2, 
there  is  a  fundamental  change  in  the  structure  of  the  OH  concentrations  with 
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Fig.  6.  Comparison  of  OH 
reaction  rate  fields  WqhU.O 
in  [mol/cm^-sl,  obtained 
from  the  conserved  scalar 
imaging  data  via  the  laminar 
flamelet  construction  proce¬ 
dure.  Shown  are  results  for 
Cases  2  (left)  and  3  (right). 
Table  gives  relative  strain 
rates  for  each  case  at  each 
downstream  location. 
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increasing  downstream  distance  from  the  nozzle.  In  particular,  the  layers 
containing  high  OH  concentrations  evident  at  the  flame  tip,  and  to  a  lesser 
extent  also  evident  at  x/L  -  0.75,  are  virtually  gone  at  x/L  -  0.50,  and  by  x/L  -  0.25 
the  OH  concentration  field  shows  only  featureless  blobs  of  high  OH  values.  Case 
3,  by  comparison,  shows  little  evidence  of  any  layer-like  structure  in  the  OH 
field  anywhere  in  the  flame.  These  fields  resemble  far  more  the  blobs  seen  at 
x/L  -  0.25  and  0.50  in  Case  2.  Indeed,  the  strain  rates  in  the  four  panels  shown 
for  Case  3  are  all  roughly  comparable  to  or  higher  than  those  at  x/L  -  0.25  in 
Case  2.  Even  at  the  flame  tip  in  Case  3,  the  strain  rates  are  still  nearly  as  high  as 
those  at  x/L  -  0.25  in  Case  2.  Similarly,  when  we  compare  panels  at  the  same  x/L 
values  between  Cases  2  and  3,  there  is  a  clear  tendency  for  the  layers  present  in 
Case  2  to  change  to  blobs  owing  to  the  increased  depth  of  nonequilibrium. 

What  is  remarkable,  however,  is  that  despite  the  wide  variations  in  the 
OH  concentration  field  structure  seen  in  Fig.  5  with  changing  downstream 
position  and  with  global  Damkohler  number,  all  of  these  fields  ultimately 
resulted  from  the  same  conserved  scalar  data  plane  shown  in  Fig.  1.  Thus  even 
the  broad  regions  of  high  OH  concentration  in  Case  3,  and  at  small  x/L  in  Case  2, 
all  resulted  from  a  simple  layer-like  structure  in  the  scalar  dissipation  rate  field  in 
Fig.  lb.  This  appears  to  have  considerable  potential  for  organizing  the  wide 
variations  apparent  in  the  OH  concentration  fields  into  a  single  underlying 
canonical  structure.  Evidently,  the  simple  layer-like  structure  in  Fig.  1  suffices  to 
reconcile  both  the  layer-like  and  broad  OH  zones  appearing  in  Fig.  5  under 
differing  conditions  of  nonequilibrium  and  locations  in  the  flame. 

Figure  6  now  shows  the  OH  reaction  rate  field  for  the  same  two  cases. 
Black  now  corresponds  to  zero  reaction  rate,  while  yellow  through  red 
correspond  to  increasing  OH  formation  rates,  and  blue  through  green  denote 
increasing  rates  of  consumption  of  OH,  as  indicated  by  each  color  bar.  Black 
regions  inside  the  flame  correspond  to  strained  out  areas.  There  are  two  main 
points  to  note  here.  First,  there  is  a  fair  correspondence  between  the  OH 
concentration  and  reaction  rate  fields  up  to  about  x/L  -  0.5,  but  beyond  that  the 
OH  mass  fraction  gives  only  a  rather  poor  indication  of  the  reaction  zone 
structure.  Second,  while  there  are  clearly  reactions  occurring  in  blobs  in  parts 
of  the  flame,  there  appears  to  be  a  persistent  layer-like  structure  in  the  reaction 
rate  field,  even  when  the  OH  concentration  field  gives  no  evidence  of  it,  as  at  x/L 
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Fig.  8.  Comparison  of  NO 
reaction  rate  fields  w^^(x,() 
in  [mol/cm^-sl,  obtained 
from  the  conserved  scalar 
imaging  data  via  the  laminar 
flamelet  construction  proce¬ 
dure.  Shown  are  results  for 
Cases  2  (left)  and  3  (right). 
Table  gives  relative  strain 
rates  for  each  case  at  each 
downstream  location. 
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-  1  in  Case  3- 

Figures  7  and  8  show  corresponding  results  for  the  NO  mass  fraction  and 
reaction  rate  fields.  (Our  work  on  nitric  oxide  in  turbulent  diffusion  flames  also 
involves  funding  from  the  Gas  Research  Institute.)  Here  it  is  necessary  to  plot 
these  in  logarithmic  form,  owing  to  the  wide  range  of  concentrations  and 
reaction  rates  present.  Similar  patterns  can  be  seen  as  the  degree  of 
nonequilibrium  increases  and  the  position  in  the  flame  changes.  In  this  case, 
however,  the  relatively  broad  region  of  ^  covered  by  the  NO  concentrations, 
coupled  with  the  strong  dependence  on  the  dissipation  rate,  leads  to  features 
not  seen  in  the  OH  fields  in  Figs.  5  and  6. 

The  flamelet  approach  demonstrated  here  for  constructing  instantaneous 
species  concentration  and  reaction  rate  fields  under  increasing  degrees  of 
nonequilibrium  chemistry  allows  considerable  insights  to  be  gained  into  the 
structure  of  turbulent  flames.  The  results  show  that  the  wide  variations  apparent 
in  the  structure  of  OH  concentration  fields,  both  in  the  results  presented  here 
and  in  direct  OH  imaging  measurements,  can  be  reconciled  into  a  single 
underlying  canonical  structure.  In  particular,  the  simple  layer-like  structure  in 
the  scalar  and  dissipation  rate  fields  seen  in  Fig.  1  suffices  to  explain  both  the 
layer-like  and  broad  OH  zones  appearing  in  Fig.  5  under  differing  conditions  of 
nonequilibrium  and  at  various  locations  in  the  flame.  Analysis  of  this  type  on  our 
high-resolution  conserved  scalar  imaging  data  are  providing  new  insights  into 
the  structure  of  chemical  reactions  in  turbulent  combustion. 

1.3.  Topology  and  Scaling  Properties  of  Scalar  Dissipation  Layers 

Molecular  mixing  of  dynamically  passive  conserved  scalar  quantities  in 
turbulent  flows  plays  a  central  role  in  an  enormous  range  of  practical  problems. 
The  mixing  process  can  be  formally  quantified  in  terms  of  the  conserved  scalar 
field  and  the  associated  instantaneous  scalar  energy  dissipation  rate  field 

{ReScY  'v;.vc(x,o.  This  aspect  of  our  work  involves  results  obtained  using  a 
laser  imaging  diagnostic  specifically  designed  for  very  highly  resolved,  four¬ 
dimensional  measurements  of  the  full  space-  and  time-varying  conserved  scalar 
field  and  the  associated  scalar  energy  dissipation  rate  field  V^  V^(x,r)  in 

turbulent  flows.  Measurements  with  this  technique  have  been  obtained  in  the 


11 


AFOSR  -  89  -  0541 

The  University  of  Michigan 


1991  -  92  Annual  Report 


self-similar  far  field  of  axisymmetric  turbulent  jets  at  outer-scale  Reynolds 
number  Re^  =  (u5/v)  as  high  as  6000,  with  the  resolution  achieved  being  finer 
than  the  local  strain-limited  molecular  diffusion  scale  A.i>  Since  the  imaged 
region  of  the  flow  is  quite  small  in  comparison  with  the  local  outer  scale  5,  and 
comparable  to  the  inner  scale  of  the  flow,  the  resulting  measured  turbulent 
scalar  field  fine  structure  should  be  essentially  independent  of  the  Reynolds 
number  and  the  flow,  and  depend  only  on  the  Schmidt  number.  In  this  sense, 
many  features  of  the  fine  structure  captured  within  our  four-dimensional  data 
space  should  to  a  large  degree  be  generic  to  large  Schmidt  number  mixing  in  ail 
high  Reynolds  number  turbulent  flows. 

The  measurements  are  of  the  aqueous  concentration  of  a  dynamically 
passive  dilute  laser  fluorescent  dye  (disodium  fluorescein)  carried  by  one  of  the 
fluids  in  the  self-similar  far  field  of  an  axisymmetric  turbulent  jet.  The  dye 
mixture  fraction  is  a  conserved  scalar  variable  for  which  Sc  =  2075.  This  mixture 
fraction  is  measured  repeatedly  in  time  throughout  a  small  three-dimensional 
volume  in  the  flow  (see  Fig.  1)  by  imaging  the  laser  induced  fluorescence  from 
dye-containing  fluid  in  the  path  of  a  laser  beam  rapidly  swept  in  a  raster  fashion 
through  the  volume  onto  a  high-speed  planar  photosensitive  array.  A  pair  of 
very  low  inertia,  galvanometric  mirror  scanners  are  used  to  synchronously 
sweep  a  collimated  laser  beam  in  a  raster  scan  fashion  through  the  desired 
volume  in  the  flow  field.  The  resulting  laser  induced  fluorescence  intensity  is 
measured  with  a  256  x  256  imaging  array,  having  center-to-center  pixel  spacings 
of  40(im.  The  array  is  synchronized  to  the  same  clock  that  drives  the  scanners, 
and  can  be  driven  at  variable  pixel  rates  up  to  11  MHz,  allowing  measurement  of 
successive  data  planes  at  a  continuous  rate  in  excess  of  140  planes  per  second. 
The  fluorescence  data  from  the  array  is  serially  acquired  through  a 
programmable  digital  port  interface,  digitized  to  8-bits  digital  depth,  then 
routed  into  a  16  MB  high-speed  dual-ported  data  buffer  from  which  it  can  be 
continuously  written  in  real  time  to  a  3-1  GB  high-speed  parallel  transfer  disk 
rank.  The  overall  sustained  data  throughput  rate  to  the  disks,  accounting  for  all 
line  and  frame  overhead  cycles,  is  up  to  9.3  MB/sec.  The  3.1  GB  disk  capacity 
can  accommodate  more  than  50,000  such  measured  256  x  256  spatial  data  planes 
within  the  four-dimensional  spatio-temporal  data  space.  The  spatial  separation 
between  adjacent  points  within  each  data  plane,  and  between  adjacent  data 
planes  within  each  data  volume,  is  smaller  than  the  local  strain-limited  molecular 
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diffusion  lengthscale  Xp  of  the  scalar  field.  Similarly,  the  temporal  separation 
between  adjacent  data  planes  within  each  data  volume,  and  between  the  same 
data  plane  in  successive  data  volumes,  is  shorter  than  the  local  molecular 
diffusion  scale  advection  time  This  resolution,  together  with  the  high  signal 
quality  attained,  allows  accurate  differentiation  of  the  measured  conserved  scalar 
data  in  all  three  space  dimensions  and  in  time  to  determine  the  components  of 
the  local  instantaneous  scalar  gradient  vector  field  V^(x,r)  at  every  point  in  the 
data  space.  From  this,  the  true  instantaneous  scalar  energy  dissipation  rate  field 
can  be  determined. 

Since  the  principal  interest  is  in  obtaining  the  scalar  energy  dissipation 
rate  field  from  the  measured  conserved  scalar  field  via  direct  differentiation  of 
the  data,  the  central  issue  is  the  spatial  and  temporal  resolution  achieved  by  the 
measurements.  From  the  measured  thickness  of  the  imaged  portion  of  the  laser 
beam,  together  with  the  pixel  size  and  the  image  ratio  of  the  measurements,  the 
volume  in  the  flow  (Ax  Ay  Az)  imaged  onto  each  pixel  can  be  determined. 
Furthermore,  for  the  pixel  clock  rates  used,  the  time  At  between  acquisition  of 
successive  data  planes  within  each  spatial  data  volume,  and  the  time  AT  between 
the  same  data  plane  in  successive  data  volumes,  can  also  be  determined.  To 
assess  the  resulting  relative  resolution  achieved,  these  smallest  spatial  and 
temporal  scales  discernible  in  the  data  must  be  compared  with  finest  local 
spatial  and  temporal  scales  on  which  gradients  in  the  conserved  scalar  field  can 
be  sustained  in  the  flow.  In  particular,  for  local  diffusion  of  vorticity  in  the 
presence  of  a  time-varying  strain  rate  £,  the  competing  effects  of  strain  and 
diffusion  establish  an  equilibrium  strain-limited  vorticity  diffusion  layer  thickness 
Xy  -  (v/e)^'2,  closely  related  to  the  Kolmogorov  scale,  giving  the  finest  scale  on 
which  spatial  gradients  in  the  strain  rate  and  vorticity  fields  can  be  locally 
sustained  in  the  flow  (e.g.  Burgers  1948;  Townsend  1951).  A  similar  competition 
between  the  effects  of  strain  and  molecular  diffusion  of  the  conserved  scalar 
establishes  a  local  strain-limited  molecular  diffusion  layer  thickness  X-d  ~(D/e)'^2^ 
related  to  the  Batchelor  scale  and  giving  the  smallest  scale  on  which  spatial 
gradients  in  the  conserved  scalar  field  can  be  sustained  by  the  flow.  The  ratio  of 
the  vorticity  and  scalar  diffusivities,  v  and  D  respectively,  establishes  the  relation 
between  these  two  scales  as  Xq  =  Note  that,  with  the  highest  strain 

rates  occurring  locally  in  the  flow  scaling  as  e  -  (u/5)-/?e*^2  the  strain-limited 
diffusion  scale  in  the  conserved  scalar  field  is  (Xj^5)  ~  The  most 
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recent  measurements  give  the  resulting  proportionality  constant  as  11.2. 

Figure  1  shows  the  conserved  scalar  data  in  a  typical  256  x  256  spatial  data 
plane.  All  quantities  presented  have  been  nondimerrsionalized  by  the  inner- 
scale  reference  values  1*  u*  =  (v/X^)  and  the  local  mean  scalar  value  s 
It  also  shows  the  true  scalar  energy  dissipation  rate  field  V^  V^(x,r)  obtained  by 
direct  differentiation  of  the  data  in  the  three  adjacent  scalar  planes  centered  on 
each  plane.  Linear  central  difference  approximations  have  been  used  to  evaluate 
the  three  components  of  the  scalar  gradient  vector  field  V^(x,r),  with  no  explicit 
smoothing  or  filtering  of  the  results.  Owing  to  the  wide  range  of  resulting 
dissipation  rates,  the  logarithmic  compression  V^-V^(x,r)  in  Fig.  1  allows  the 
structure  at  low  dissipation  rates  to  be  examined.  Our  previous  work  on  these 
fine  structure  maps  of  the  instantaneous  scalar  dissipation  field  has  clearly  shown 
that  essentially  all  of  the  molecular  mixing  occurs  in  thin  sheet-like  layers  of  the 
type  described  by  Burgers  and  Townsend  for  the  vorticity  field.  The  internal 
structure  of  the  molecular  mixing  within  the  layers  apparent  in  these 
measurements  can  also  be  examined  from  our  measurements,  and  confirm  that 
strain-limited  solutions  of  the  Burgers  and  Townsend  form  give  a  remarkably 
accurate  description  of  the  true  scalar  energy  dissipation  profiles  within  these 
layers. 


Our  work  with  these  measurements  and  the  resulting  data  during  the  past 
twelve  months  has  focused  on  analyses  using  three  different  analytical 
approaches.  These  include  (1)  layer  separation  distributions,  (2)  critical  point 
analyses,  and  (3)  multifractal  topology  studies.  Our  interests  are  principally  in 
the  topology  which  the  scalar  dissipation  layers  assume  as  a  result  of  their 
repeated  stretching  and  folding  by  the  turbulent  strain  rate  and  vorticity  fields. 
The  aim  is  to  obtain  measurements  of  the  scalings  which  this  topology  satisfies, 
and  relate  these  to  the  dynamics  of  the  “turbulent  cascade”  represented  by  this 
stretching  and  folding  process.  We  are  motivated  by  recent  successes  in 
comparatively  simple  two-dimensional  chaotic  flows  (Ottino  1992;  Muzzio, 
Meneveau,  Swanson  &  Ottino  1992),  where  a  periodic  stretching  and  folding  of 
the  scalar  layers  leads  to  scalings  for  the  distribution  of  dissipation  layer 
separations  that  can  be  reconciled  with  simple  physical  arguments.  These 
analyses  suggest  a  lognormal  distribution  for  the  distribution  of  layer 
separations,  though  on  quite  different  (and  more  rigorous)  grounds  than  the 
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(X/X.) 


Fig.  9.  The  measured  distribution  of  the  scalar  dissipation  layer  separation  dis¬ 
tances  obtained  from  three-dimensional  spatial  data  volumes  of  the  type 
shown  in  Fig.  1.  Note  the  quasi-lognormal  form  obtained  for  small  layer  separa¬ 
tion  values,  and  the  non-lognormal  form  found  at  larger  separation  values,  (a)  The 
linear  version,  (b)  The  logarithmic  version. 
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classical  lognormal  predictions  for  other  single-sided  quantities  in  turbulence. 
From  the  three-dimensional  spatial  character  of  the  turbulent  dissipation  rate 
field  in  any  given  data  volume,  we  constmct  the  surface  of  local  layer-normal 
dissipation  maxima.  We  can  then  compute  the  distribution  of  dissipation  layer 
separations  The  roughly  lognormal  form  appears  consistent  with  the 

scalings  obtained  in  the  simple  periodic  flows,  but  for  large  layer  separations  we 
find  deviations  from  this  simple  scaling.  Our  analyses  of  these  differences  is  not 
yet  complete,  but  they  appear  to  be  due  to  differences  between  turbulence  in 
closed  and  open  systems.  In  particular,  in  an  open  system  like  the  turbulent  jet 
(and  other  practical  free  shear  flows),  there  is  a  continual  entrainment  of 
dissipation-free  fluid  into  the  flow.  This  dissipation-free  fluid  is  subsequently 
intertwined  with  the  dissipation  layers  in  the  turbulent  flow,  and  continually 
introduces  regions  with  large  dissipation  layer  separation  distances.  No  such 
mecharrism  exists  in  a  closed  turbulent  flow,  since  there  is  no  entrainment  of 
dissipation-free  fluid.  As  a  consequence,  the  lognormal  distribution  (which  does 
not  take  into  account  such  entrainment)  produces  a  relatively  good  fit  to  the 
distribution  for  small  separations.  For  large  values  it  would  not  be  expected  to 
give  accurate  predictions,  and  that  is  indeed  seen  to  be  the  case.  At  these  larger 
values,  the  scaling  appears  to  be  consistent  with  simple  volume-preserving 
constraints  that  play  a  dominant  role  in  the  entrainment  process. 

These  analyses  play  a  key  role  in  developing  a  practical  model  for  reacting 
turbulent  flows  to  aid  in  development  of  new  airbreathing  propulsion  systems. 
In  particular,  the  two-parameter  nonequilibrium  reaction  chemistry  formulation 
in  §1.2  requires  joint  information  about  the  conserved  scalar  and  its  dissipation 
rate.  When  that  formulation  is  coupled  with  these  simple  scalings  for  the 
dissipation  layer  separations,  it  appears  to  be  possible  to  predict  the  relative 
dominance  of  layer-like  reaction  zone  structure  versus  broad  distributed 
reaction  zone  structure.  Moreover,  it  appears  possible  to  fold  this  information 
into  a  new  class  of  sub-grid  scale  models  for  reacting  flows,  in  which  the 
computed  scales  impose  the  local  outer  scale  strain  field  on  the  distribution  of 
layer  separations.  This  scales  the  distribution  of  local  strain  rates  at  the  subgrid 
scale,  which  in  turn  coupled  with  the  nonequilibrium  chemistry  formulation  to 
allow  prediction  of  species  mass  fractions  and  reaction  rates  for  complex 
reaction  chemistry. 
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2.  Personnel 

The  following  personnel  were  associated  with  this  research  effort  during  the 
report  period. 

A.  Faoilty 

Dr.  Werner  J.A.  Dahm;  Associate  Professor,  The  University  of  Michigan, 
Department  of  Aerospace  Engineering,  Ph.D.,  US  Citizen. 

B.  Graduate  Students 

Southerland,  K.B.;  Graduate  Student,  The  University  of  Michigan,  Department  of 
Aerospace  Engineering,  Ph.D.  program,  US  Citizen.  Mr.  Southerland  is  the 
principal  graduate  student  research  assistant  working  on  this  grant,  and  he  is  fully 
supported  through  his  work  on  the  grant. 

Bish,  E.;  The  University  of  Michigan,  Department  of  Aerospace  Engineering, 
Ph.D.  program,  US  Citizen.  Mr.  Bish  is  85%  supported  through  a  DoEd 
Fellowship;  he  receives  15%  support  ($200/month  +  1  summer  month)  from  his 
work  on  this  grant. 

Frederiken,  R.E.;  The  University  of  Michigan,  Department  of  Aerospace 
Engineering,  Ph.D.  program,  US  Citizen.  Mr.  Frederiken  is  80%  supported 
through  a  Departmental  Teaching  Assistantship;  he  receives  10%  support 
($100/month  +  1.5  summer  months)  from  his  work  on  this  grant.  Application  to 
the  DoD  ASSERT  Program  has  been  made  to  allow  him  to  work  full-time  on  this 
project  with  funds  outside  the  grant  budget. 

C.  Post-Doctorates 

Buch,  Kenneth  A.,  Jr.;  Michigan,  Aerospace  Engineering,  Ph.D.,  US  Citizen. 
During  the  time  spanned  by  this  worksheet.  Dr.  Buch  worked  partly  on  this  grant 
from  June  1991  through  December  1991.  He  received  approximately  20%  of  his 
support  f  his  work  on  this  grant. 
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3.  Publications 


The  following  publications  have  resulted  during  this  reporting  period  in 
part  from  the  research  work  undertaken  in  this  research  program  into  the 
structure  of  chemical  reactions  in  turbulent  flows. 

1.  Dahm,  W.J.A.  &  Bish.  E.  (1992)  “High  resolution  measurements  of  molecular 
transport  processes  in  turbulent  combustion,”  Turbulence  and  Molecular  Proce.sses 
in  Combustion.  Elsevier,  in  press. 

2.  Dahm,  W.J.A.,  Southerland,  K.B.  and  Buch,  K.A.  (1991)  Four-dimensional  laser 
induced  fluorescence  measurements  of  conserved  scalar  mixing  in  turbulent  flows. 
Applications  of  Laser  Techniques  to  Ruid  Mechanics.  3-18,  R.  Adrian,  Ed., 
Springer  Verlag,  Berlin. 

3.  Dahm,  W.J.A.  &  Southerland,  K.B.  (1992)  “Quantitative,  four-dimensional,  laser 
induced  fluorescence  imaging  studies  of  scalar  mixing  in  turbulent  flows,”  Flow 
Visualization  6. 243-247,  Springer- Verlag,  Berlin. 

4.  Buch,  K.A.,  Dahm,  W.J.A.,  Dibble,  R.W.  &  Barlow,  R.S.  (1992)  “Equilibrium 
structure  of  reaction  rate  fields  in  turbulent  jet  diffusion  flames,”  to  appear  in 
Proceedings  of  the  24th  International  Symposium  on  Combustion.  The  Combustion 
Institute,  Pittsburgh. 

5.  Dahm,  W.J.A.  (1992)  “Experimental  studies  of  the  fine  scale  structure  of  mixing  in 
turbulent  flows,”  Proceedings  of  the  13th  Symposium  on  Turbulence.  The 
University  of  Missouri  (Rolla),  Rolia,  MO. 


4.  Presentations 


The  following  presentations  have  resulted  during  this  reporting  period  in 
part  from  the  research  work  undertaken  in  this  research  program  into  the 
structure  of  chemical  reactions  in  turbulent  flows. 


1.  Fully-Resolved,  Four-Dimensional,  Laboratory  Measurements  of  the  Fine  Scale 
Structure  of  Turbulent  Flows  -  Invited  Seminar,  Mechanical  Engineering 
Department,  Purdue  University,  West  Lafayette,  IN,  September  1991. 

2.  Fully-Resolved,  Four-Dimensional,  Laboratory  Measurements  of  the  Fine  Scale 
Stmcture  of  Mixing  in  Turbulent  Flows  -  Invited  Seminar,  Mechanical  Engineering 
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Department,  University  of  California  at  Davis,  Davis,  CA,  February  1992. 

3.  Experimental  Studies  of  the  Fine  Structure  of  Mixing  in  Turbulent  Flows  -  Invited 
Seminar,  Mechanical,  Aerospace  &  Nuclear  Engineering  Department,  University  of 
California  (Los  Angeles),  March  1992. 

4.  Small  Scale  Structure  of  Turbulent  Diffusion  Flames  -  Invited  Lecture,  19th  Meeting 
of  the  Sandia  Technical  Group  on  Aerothermochemistry  of  Reacting  Flows, 
Lawrence  Berkeley  Laboratories,  CA,  April  1992. 

5.  Structure  of  Mixing  and  Reactions  in  Turbulent  Rows  -  Invited  Seminar,  High 
Temperature  Gas  Dynamics  Laboratories,  Mechanical  Engineering  Department, 
Stanford  University,  Stanford,  CA,  April  1992. 

6.  Equilibrium  Structure  of  Chemical  Reaction  Rate  Fields  in  Turbulent  Diffusion 
Rames  -  24th  International  Symposium  on  Combustion,  The  University  of  Sydney, 
Sydney,  Australia,  June  1992. 

7.  Experimental  Studies  of  the  Fine  Scale  Structure  of  Mixing  in  Turbulent  Rows  - 
Invited  Lecture,  Thirteenth  Symposium  on  Turbulence,  Rolla,  MO,  September 
1992. 

8.  Quantitative,  Four-Dimensional,  Laser  Induced  Ruorescence  Imaging  Studies  of 
Scalar  Mixing  in  Turbulent  Rows  -  Sixth  International  Symposium  on  Row 
Visualization,  Yokohama,  Japan,  October  1992. 

9.  High  Resolution  Measurements  of  Molecular  Transport  Processes  in  Turbulent 
Combustion  -  Invited  Lecture,  6th  Toyota  Symposium  on  Turbulence  and  Molecular 
Processes  in  Combustion,  Nagoya,  Japan,  (October  1992. 

10.  Modeling  NOx  in  Turbulent  Diffusion  Rames  -  Invited  Presentation,  2nd 
International  Workshop  on  NOx  in  Natural  Gas  Flames,  Orlando,  FL,  November 
1992. 

1 1 .  Laboratory  Studies  of  the  Fine  Scale  Structure  of  Turbulent  Rows:  5c  »  1,  5c  «  1 , 
and  Implications  for  Reacting  Rows  -  Invited  Seminar,  Theoretical  Division  TA-3, 
Los  Alamos  National  Laboratories,  December  1992. 


5.  Research-Related  Interactions 

Interactions  related  to  this  research  effort  during  the  reporting  period 
have  been  underway  with  personnel  at  the  Combustion  Research  Facility  (CRF) 
at  Sandia  National  Laboratories  (Livermore,  CA);  at  ETH-Zurich  in  Switzerland; 
and  at  Nagoya  University  in  Japan.  The  Principal  Investigator  conducted  a  five- 
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month  sabbatical  visit  to  the  CRF  from  January  -  May  1992  during  which 
measurements  of  OH  concentrations  were  made  in  conjunction  with  Dr.  P.H. 
Paul  for  comparison  with  our  flamelet  analysis  predictions.  Measurements  were 
also  planned  for  four-dimensional  Rayleigh  imaging  of  a  conserved  scalar  field  in 
high  Reynolds  number  turbulent  mixing.  An  on-going  interaction  with  the  Swiss 
Federal  Institute  of  Technology  (ETH)  in  Zurich  has  been  aimed  at  setting  up  a 
similar  diagnostic  system  as  ours  in  their  laboratories  for  additional 
measurements.  A  graduate  from  our  research  group  has  been  sent  to  ETH  on  a 
postdoctoral  assignment  for  this  purpose.  A  visit  to  our  laboratories  by 
Professor  Franz  Rys  of  ETN-Lausanne  was  made  as  part  of  this  collaboration.  An 
ongoing  interaction  with  Professor  T.  Takeno  of  Nagoya  University  is  aimed 
largely  at  further  investigations  of  flamelet  chemistry. 


6.  Summary  of  Inventions 

No  invention  disclosures  are  to  be  filed  on  the  basis  of  work  done  during 
the  reporting  period.  The  invention  disclosure  filed  during  the  previous 
reporting  period  is  still  being  evaluated  by  the  University’s  patent  counsel. 
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